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Abstract 

We present a weak-lensing analysis for the merging galaxy cluster, CIZA J2242.8+5301, hosting double 
radio relics, using three-band Subaru/Suprime-Cam imaging {Br'z'). Since the lifetime of dark matter halos 
colliding into clusters is longer than that of X-ray emitting gas halos, weak-lensing analysis is a powerful 
method to constrain a merger dynamics. Two-dimensional shear fitting using a clean background catalog 
suggests that the cluster undergoes a merger with a mass ratio of about 2:1. The main halo is located 
around the gas core in the southern region, while no concentrated gas core is associated with the northern 
sub halo. We find that the projected cluster mass distribution resulting from an unequal-mass merger is 
in excellent agreement with the curved shapes of the two radio relics and the overall X-ray morphology 
except for the lack of the northern gas core. The lack of a prominent radio halo enables us to constrain 
an upper limit of the fractional energy of magneto-hydrodynamics turbulence of {5B/BY < 0(10“®) at a 
resonant wavenumber, by a balance between the acceleration time and the time after the core passage or 
the cooling time, with an assumption of resonant acceleration by second-order Fermi process. 

Keywords: galaxies: clusters: individual (CIZA J2242.8+5301) - gravitational leasing: weak-X-rays: 
galaxies: clusters 


1. Introduction 

The acceleration mechanism generating relativistic par¬ 
ticles in the intracluster medium (ICM) along with its 
associated non-thermal phenomena is one of the long¬ 
standing unresolved problems in studies of cluster physics. 
The presence of relativistic electrons and magnetic fields 
is recognized by diffuse, synchrotron radio emission (on 
Mpc scales), which is not associated with compact sources 
such as AGN, but with cluster merging (e.g. Feretti et al. 
2012; Brunetti & Jones 2014). The diffuse radio emission 
in galaxy clusters is classified according to their morphol¬ 
ogy into two types, radio halos and relics. Radio halos, to 
first order, follow the X-ray brightness distribution in the 
central region, sometimes along with a curved spectrum 
(e.g. Brunetti et al. 2008). Radio relics display a filamen¬ 
tary structure in the peripheral region. Some relics have 
a harder spectrum than that of radio halos (e.g. Brunetti 
et al. 2008; van Weeren et al. 2010). The fundamental 

* Based on data collected at Subaru Telescope, which is operated 
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physics to form radio halos and relics and give rise to dif¬ 
ferent morphological types are one of the most outstand¬ 
ing problems. Furthermore, the acceleration mechanism is 
still unclear from the theoretical side. In theory, two types 
of acceleration mechanism have been proposed, the first 
(e.g. Bell 1978; Blandford & Ostriker 1978; Drury 1983) 
and second order Fermi acceleration (e.g. Schlickeiser et al. 
1987; Ohno et al. 2002; Fujita et al. 2003; Brunetti et al. 
2004). It is unresolved through what process relativistic 
particles are accelerated/re-accelerated in merging clus¬ 
ters. Since the merger dynamics is controlled by dark 
matter, it is of vital importance for understanding non- 
thermal phenomena to measure the cluster mass distribu¬ 
tion for merging clusters, (e.g. Okabe & Umetsu 2008; 
Okabe et al. 2010, 2011; Medezinski et al. 2015). 

Recent radio observation (van Weeren et al. 2010) dis¬ 
covered spectacular, double-giant radio relics in CIZA 
J2242.8+5301. The northern relic is larger than the 
southern one, which suggests non-equal-mass merger (van 
Weeren et al. 2011). The northern giant relic shows a 
spatial distribution of aging relativistic electrons, indicat¬ 
ing a presence of strong magnetic field (5 —7^G). The 
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magnetic field at the radio relic is stronger than equipar- 
tition magnetic fields in radio halos (~ 0.1 — 1 /rG), but 
consistent with a lower limit of the magnetic field (> 3/iG) 
in other relics by comparing the inverse Compton X- 
ray emission with the measured radio synchrotron emis¬ 
sion (e.g. Nakazawa et al. 2009; Finoguenov et al. 2010). 
Assuming a diffusive shock acceleration (DSA; Drury 
1983), the Mach number of shock velocity is expected to 
he M = 4. 6 ^ 0 , 9 - Suzaku X-ray observation (Akamatsu 
& Kawahara 2013) has discovered a temperature jump 
across the northern relic and estimated that the Mach 
number is somewhat lower than expected by radio obser¬ 
vation. A discrepancy between the X-ray-estimated and 
radio-estimated Mach numbers makes it puzzling to un¬ 
derstand a particle acceleration process. On the other 
hand, Stroe et al. (2014) have pointed out that an intrinsic 
beam-image with of the radio telescope affects the estima¬ 
tion of the radio spectral index and revised Jhi = 2.901q 
which is close to X-ray measurements. Therefore, the 
study on a comparison of Mach numbers is still contro¬ 
versial. No prominent radio halo in the central region is 
found by multi-bands radio analysis (van Weeren et al. 
2010; Stroe et al. 2013). Recent study of weak-lensing 
analysis (Jee et al. 2015) has shown a bimodal structure 
of the dark matter and member galaxies and concluded 
an equal-mass merger. 

We here report weak-lensing analysis using data from 
Subaru/Suprime-Gam. The paper is organized as fol¬ 
lows: Section 2 presents the weak-lensing analysis of 
Subaru/Suprime-Gam. We present distribution of mem¬ 
ber galaxies and mass measurements in Sections 3 and 
4, respectively. In Section 5, we compare the projected 
mass distribution with the X-ray emitting gas, radio relics 
and member galaxy distribution. The results are dis¬ 
cussed in Section 6 and summarized in Section 7. We 
use = 0.3, Da = 0.7 and Hq = lOO/i km 
{Ho = 70/170 km s“^Mpc“^). 

2. Data Analysis 

We took three-band imaging data {B, r', z') for 
CIZAJ2242.8-1-5301 {zi = 0.192) with the Suprime-cam 
(Miyazaki et al. 2002) at the Subaru 8.2-m telescope. 
Here, the subscript I denotes the lensing (cluster) red- 
shift. The exposure times for the B- r'- and z'- bands 
are 24, 36 and 12 mins, respectively. Those imaging are 
taken through the Service program (P.I.:N. Okabe). The 
r'- image, of which seeing is 0.7 arcsec, is used to measure 
ellipticities of galaxies. The other bands are combined 
with the r'-band data to minimize a contamination of un- 
lensed member galaxies in background source catalog. 

We use the standard pipeline reduction software for the 
Suprime-Cam, SDFRED (Yagi et al. 2002; Ouchi et al. 
2004) modified for the new CCDs. The data reduction 
follows Okabe et al. (2013, 2014a). The zero point of mag¬ 
nitude is corrected for the Galactic extinction (Schlafly & 
Finkbeiner 2011) using a single function over the whole 
field and calibrated by comparing magnitudes of stellar 
objects with those of standard stars 
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We measure ellipticity of each galaxy in terms of weight¬ 
ing the surface brightness, following the KSB method 
(Kaiser et al. 1995) with some modifications (Okabe et al. 
2014a). The details are described in Okabe et al. (2014a). 
The photometric redshift of each galaxy is estimated as 
the ensemble average of the COSMOS photometric red- 
shift (Ilbert et al. 2013) of neighboring COSMOS galaxies 
in the three-magnitudes plane. We adopt 100 COSMOS 
galaxies for the estimation and confirm that the change by 
changing the sampling number is negligible. The galaxies 
for shape measurements are selected by rh > r'^ + Srh*, 
fg > fg + Sr* and 22 < r', where and rg are the half- 
light radius and the Gaussian radius, respectively. The 
asterisk denotes the stellar objects. Since the cluster is 
located at the Galactic plane, the number density of back¬ 
ground galaxies is small 9 arcmin“^ because of a consider¬ 
able number of contaminated stars, compared to those in 
previous studies (e.g. Okabe & Umetsu 2008; Okabe et al. 
2010, 2013; Umetsu et al. 2014). 

A secure background selection is critically important 
to measure accurately cluster and subcluster masses (e.g. 
Broadhurst et al. 2005; Okabe et al. 2013; Okabe & Smith 
2015). We use color information [r' — z' and B — r') for the 
selection (e.g. Medezinski et al. 2010; Umetsu et al. 2010). 
The red-sequence of member E/SO galaxies appeared in [r' 
and r' — z') and (r' and B — r') planes is fitted with a linear 
function. The mean lensing depth (Dis/Ds) is computed 
as a function of the color offsets from the red-sequence 
of (U - z') - (r' - z')e/so and {B - r') - {B - r')E/so- 
Here, and Dis are the angular diameter distances 
to the sources and between the cluster (lens) and the 
sources, respectively. The resulting (Dis/Dg) map is 
shown in the left panel of Eigure 1. The grid box size 
is O.lABmag. (Dis/Dg) at small color-offsets presents low 
lensing depth, suggesting a finite contamination of mem¬ 
ber galaxies. Eollowing Medezinski et al. (2010), we first 
select galaxies of which {B — r') is bluer than those of 
clusters {{B — r') — {B — r')E/so < —o'B-r'), where as-r' 
is the color width of the red sequence. In a similar manner, 
galaxies of which colors are located in the width of cr^'-z' 
are excluded in our source background catalog. Since a 
combination of the lensing depth and the color informa¬ 
tion is a good estimator to monitor the color distribu¬ 
tion for possible member galaxies, we exclude color regions 
for our background sources by employing {Dig/Ds) < 0.6, 
which is equivalent to (^s) < 0.52. The number density 
of background source galaxies is 2.6 arcmin“^. The mean 
source redshift is {zs) = 0.639. 

3. Galaxy Distribution 

We make a map of member galaxies for visual pur¬ 
pose to depict the configuration of cluster merger. We 
here define bright member galaxies with r' < 21ABmag, 
\{B - r') -{B - r')E/sol < and |(r' - z') - {r' - 

z')e/soI < o'r'-z'- The luminosity in the r'—band is com¬ 
puted with K-correction assuming a single cluster redshift. 
We make a luminosity map convolved with a Gaussian 
smoothing kernel oc exp[—0^/0^], where 6g = 2' corre- 
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Fig. 1. Left : (Dis/Ds) distribution in the color-offset plane of (r' — z') — {r' — ^Oe/SO (-^ ~ — ?"0e/SO- Right: black 

points are all galaxies for which shapes are measured. Blue points are background galaxies. Green region denotes the color area 
with (Dis/Ds) < 0.6. Red dashed lines are the width of the red sequence. 


spending to FWHM= 3f3. The right panel of Figure 2 
clearly shows a bimodal structure of member galaxies. 
The luminosity peaks are associated with the brightest 
cluster galaxies for two halo components, which provides 
us with priors of the number of halos and their central 
positions. 

4. Mass Measurement 

We simultaneously fit a two-dimensional shear pattern 
with multi-components model, following previous papers 
(see details; Okabe et al. 2011; Watanabe et al. 2011; 
Merluzzi et al. 2015). We pixelize the shear pattern into a 
regular grid of 1' x 1' without any spatial smoothing pro¬ 
cedures. We employ the NFW mass model predicted by 
numerical simulations (Navarro et al. 1996) as the mass 
model. The NFW mass density prohle is specified by 
two parameters of the mass Ma and the halo concen¬ 
tration CA, where A denotes A times the critical mass 
density, pcr(+), at the cluster redshift. We also treat 
centers (ac,Sc) of the two mass components as free pa¬ 
rameters (e.g. Oguri et al. 2010). Based on a prior of 
member galaxy distribution (Section 3; right panel of 
Figure 2), we consider two-halo components in the north¬ 
ern and southern region. There are eight parameters in 
total p = ,S^). Here, the su¬ 

perscripts denote the southern halo (S) and the north¬ 
ern halo (N), respectively. We adopt the Markov Chain 
Monte Carlo method with standard Metropolis-Hastings 
sampling. Since it is difficult to determine the concen¬ 
tration parameter of each individual halo because of the 
small number of background galaxies and the degener¬ 
acy between model parameters, we use a Gaussian prior 
that follows the halo concentration-mass (c-M) relation 


of Bhattacharya et al. (2013) derived from N-body cos¬ 
mological simulations. The Bhattacharya et al. (2013) 
c-M relation is in agreement with recent stacked leas¬ 
ing studies (Okabe et al. 2013; Okabe & Smith 2015; 
Umetsu et al. 2014, 2015). We use a flat prior for other 
parameters. The ranges of virial masses are restricted 
to 0 < Ma < 30 X 10^^ M^Mq. The centers are allowed 
within ±2f around the candidates of most luminous clus¬ 
ter galaxies. The best-fit parameters are shown in Table 
1. The best-fit mass for the southern halo is twice higher 
than that for the northern halo, though an equal-mass 
ratio cannot be ruled out within errors. The offsets be¬ 
tween lensing centers and most luminous galaxies on the 
sky plane for the southern and the northern halos are 
36^36® h“^kpc and 151^67^ h“^kpc, respectively. Thus, 
the offset is not significant for the southern main halo and 
at ^ 2a level for the northern subhalo. 

The tangential shear profiles of two mass components 
as a function of the projected distance from the best-fit 
center of the southern halo are shown in Figure 3. The 
center is set to be the best-fit value. The lensing signal at 
r ~ 8f is depressed by the northern mass halo. The best- 
fit lensing profile for the southern halo (green dashed line) 
monotonically decreases, while the best-fit lensing profile 
for the northern halo (blue dotted line) rapidly increases 
up at r ^ 8'. due to an off-centering effect (e.g. Yang et al. 
2006). The total mass model (red solid line) well describes 
the observed lensing profile. 
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Fig. 2. Right: the best-fit mass distribution computed from a simultaneous fit of two-dimensional shear pattern (Sec 4). Contour 
(red) levels are in units of the convergence, k, starting from k = 0.1 with steps of k = 0.1. Overlaid contours (white) are the radio 
emission (Stroe et al. 2013) at 1.2 GHz in units of signal-to-noise ratio [4,8,16, 32](7 rmS' Middle: X-ray image of Suzaku satellite 
in 0.5 —2.0keV (Akamatsu & Kawahara 2013; Akamatsu et al. 2015). Right : Luminosity map for bright member galaxies (Sec 3) 
computed with Gaussian smoothing FWHM= 3f3. 

Table 1. Mass estimates for the southern and northern mass components (10^^ h~^MQ) and their centroids in the units of degree, 
respectively. 
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5. A comparison between Projected 
Distributions of Radio Relic, ICM, and 
Mass 

The acceleration process of relativistic electrons emit¬ 
ting giant relics are believed to be associated with on¬ 
going merger, but the details remain unknown. Since the 
dark matter controls the merger dynamics, a comparison 
of projected distributions between the dark matter, the 
gas and the radio relics would provide us with an im¬ 
portant information for understanding the non-thermal 
physics (e.g. Okabe et al. 2010). Since the number density 
of background galaxies is small, we compute the projected 
mass distribution calculated from the best-fit parameters 
without any smoothing kernels, rather than weak-mass 
reconstruction (Okabe & Umetsu 2008) in order to avoid 
noisy feature. The best-fit projected mass distribution 
makes it clear to understand the configuration of cluster 
merger. We also investigate the projected mass distribu¬ 
tion by changing parameters with Icr measurement errors, 
the overall feature is not significantly changed. The resul¬ 
tant maps are shown in the left panel of Figure 2. The left 
panel is the best-fit mass distribution, overlaid with radio 
contours (Stroe et al. 2013). The more and less massive 
clusters are located at the south and the north, respec¬ 


tively. Both the northern and the southern relics are very 
similar to the projected mass distribution. In particular, 
the curved shapes of two relics are tightly correlated with 
the curvature of the iso-contours of the projected mass 
distribution. It indicates that the bimodal mass distribu¬ 
tion with a different mass ratio has some influence on the 
formation of asymmetric radio relics. 

The middle panel compares the X-ray surface bright¬ 
ness distribution with the mass distribution and the ra¬ 
dio relics. The X-ray distribution is elongated along 
the north-south direction. The overall gas distribution 
is in an excellent agreement with the mass distribution. 
Especially, the faint X-ray emission (green color) at ^ 5' 
from the merger axis follows the curvature of dark matter 
distribution. The X-ray peak is associated with the south¬ 
ern massive cluster, while no prominent gas core is found 
around the center of the northern less massive cluster. 
It suggests that the gas core associated with the north¬ 
ern subcluster is disrupted by the merger (Tormen et al. 
2004). 

The smoothed luminosity map for member galaxies 
(right panel; Sec. 3) shows a clear bimodal spatial dis¬ 
tribution which is similar to weak-lensing mass distri¬ 
bution. The luminosity around each central region is 
comparable to each other. We compute the mass-to- 
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Fig. 3. Top panel : tangential shear profile as a function 
of the projected cluster-centric radius from the center of the 
southern main halo determined by two-dimensional shear fit¬ 
ting. A bump in the profile signals is found at r^8'. The 
profile is well described by two NFW components of the north¬ 
ern and southern halos. The red solid, green dashed and blue 
dotted lines are the best-fitting NFW profile for the two ha¬ 
los, the southern halo and the northern halo, respectively. 
Bottom panel: the 45° rotation components, 0 -gx- 

light ratio within < 250/1^0^ with central positions de¬ 
termined by two-dimensional shear pattern. We obtain 
M/L{< 250/ifQ^) = 687 ± 245 /i7oMq/Lq for the southern 
halo and A29 ±18AhyQMQ/L q for the northern halo, re¬ 
spectively. Here, we use the projected mass for each 
halo component and the r-band total luminosity corrected 
with a Schechter luminosity function in the central region 
(< SOO/ifo^; Popesso et al. 2005). The mass-to-light ratio is 
somewhat larger than those in other clusters (e.g. Okabe 
& Umetsu 2008; Medezinski et al. 2010). This could be 
caused by the fact that some member galaxies are missing 
due to a contamination of a considerably number of stars 
in Galactic-plane position. Although the viral mass for 
the southern halo is about two times higher than that for 
the northern halo, the mass-to-light ratios are comparable 
because two masses within small radii are comparable due 
to a difference of the concentration parameter (Table 1). 

6. Discussion 

The best-lit parameters obtained by two-dimensional 
shear fitting suggest that CIZAJ2242.8-1-5301 undergoes 
the cluster merger with ~ 2 :1 mass ratio. The main clus¬ 
ter is the southern halo and the subcluster is the north¬ 
ern halo, respectively. The overall distribution of X-ray 
emitting gas is similar to the projected mass distribution 
(the middle panel in Figure 2). On the other hand, the 
gas core is associated with the main cluster, while no sig¬ 


nificant gas core around the northern subhalo is found. 
It indicates that the gas core associated with the main 
cluster survives under the collisional matter interaction 
and the gas core with the subhalo is disrupted (Tormen 
et al. 2004). Chandra X-ray observation (Ogrean et al. 
2014) has reported possible density jumps in the surface 
brightness profile at inner radii between two relics. Visual 
inspection of their result suggests that those jumps fol¬ 
low the overall mass distribution. The Chandra X-ray 
image also shows that the X-ray emission from the north¬ 
ern central region (at r 2' from from the approximate 
northern subhalo center) is somewhat more extended than 
that from the southern central region. It would point out 
that the disrupted gas core in the subhalo is more ex¬ 
tended. Although the gas distribution does not always 
follow the dark matter distribution in violent merging 
clusters (Okabe & Umetsu 2008), the gas distribution is 
similar to the dark matter distribution in the cluster on 
a very early stage after the core passage. The big giant 
radio relic is in the front of an advancing direction of the 
northern subhalo. Asymmetric radio relics are very simi¬ 
lar to the projected mass distribution. 

Since the dark matter subhalo is less disrupted by the 
cluster merger compared to the collisional gas core, the 
positions of dark matter halos provide us with the fun¬ 
damental information of the cluster merger. The dis¬ 
tance between weak-lensing determined centers of two ha¬ 
los is 712 /i“^kpc. The center of mass can be esti¬ 
mated to be (340.667deg, 53.0202deg) assuming a point 
mass. The distances between outer edge of the north¬ 
ern relic and the center of mass, c?nr-cm, and between 
outer edge of the southern relic and the center of mass, 
c^SR-CM are ~ 91 and ~ 514, respectively. Interestingly, 
the distance ratio is approximately equal to the ratio 
of the Mach numbers estimated by X-ray temperature, 
c1nr-cm : dsR-CM ^ A4nr : A4sr (Akamatsu et al. 2015), 
where A^nr = 2.7lo 4 and A4 sr = 1.7lg;3. Furthermore, 
the projected distances between the dark matter halos and 
the center of the mass, cInm-cm and dsM-CM) follow this 
ratio cInm-cm : c1sm-cm A4nr : A4sr. Since the dis¬ 
tance is equivalent to the merger timescale multiplied by 
the velocity, those positions relative to the center of mass 
are simply explained by the timescale after the core pas¬ 
sage. In other words, the formation of radio relic is possi¬ 
bly triggered after the core passage. Our mass measure¬ 
ment could consistently explain the merger configuration 
observed by X-ray, radio and weak-lensing techniques. 

The particle acceleration mechanism of relativistic elec¬ 
trons is one of the unsettled problems. A visual in¬ 
spection of the sharp morphology of radio relics (e.g. 
Brunetti & Jones 2014) and the spectral aging (van 
Weeren et al. 2010) suggest that radio relics are con¬ 
nected with merger shocks. Indeed, X-ray observations 
(Akamatsu & Kawahara 2013; Akamatsu et al. 2015) dis¬ 
covered the temperature jump around the radio relics, 
though the Mach numbers estimated by X-ray observa¬ 
tion (A4nr = 2.7 ± 0.5 and A4sr = l-f ± 0.3). are some¬ 
what lower than those (A4nr = 4.58 ± 1.09 and AdsR = 
2.81 ±0.19) by radio spectral index based on DSA model 
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(Stroe et al. 2013). Stroe et al. (2014) have revised 
7 Wnr = 2.90^;;;}° by considering a convolution of the beam 
size of radio telescopes. Weak-lensing analysis has shown 
that the morphology of radio relics is similar to the pro¬ 
jected mass distribution. It indicates that the shapes of 
radio relics on the cluster scale radio relics well follow the 
global structure of the mass distribution. The observa¬ 
tional fact implies the connection between the gas physics 
and the dark matter distribution. Two scenarios of the 
formation of radio relic are proposed. First, van Weeren 
et al. ( 2010 ) have suggested that a cluster merger shock 
accelerates electrons through some processes like the DSA 
(Drury 1983). Second, the fossil radio plasma seeded in 
the past (e.g. radio lobes from active galactic nuclei) is re¬ 
accelerated by gas compression or shocks (e.g. Enfilin & 
Gopal-Krishna 2001; Hoeft et al. 2004; Markevitch et al. 
2005). Hoeft et al. (2004) have shown by numerical sim¬ 
ulations that the radio ghost is flared at curved merger 
shocks. Although this study cannot constrain the forma¬ 
tion process of radio relics, it is worth mentioning that the 
second scenario should be discussed in the cluster com¬ 
munity, because the first scenario conflicts with other ob¬ 
servational evidences but the second does not. Prominent 
detached shocks in A520 and the bullet cluster which have 
been first discovered by X-ray observations are not associ¬ 
ated with radio relics but with radio halos (Govoni et al. 
2004). If all merger shocks injected relativistic electrons, 
the radio relic should be observed in other famous merg¬ 
ing clusters. If the formation of radio relic depends on 
the initial fossil radio plasma, we do not face the above 
problem. In order to understand the formation of radio 
relics, further systematic multi-wavelength studies of X- 
ray, radio, optical and weak-lensing analyses are vitally 
important. 

We compare with previous numerical simulations re¬ 
garding cluster mergers, van Weeren et al. (2011) have 
conducted numerical simulations for the cluster and shown 
that the northern massive cluster with the largest X-ray 
core radius triggers the largest shock wave after the core 
passage. Their result conflicts with both X-ray and weak- 
lensing measurements, because the X-ray gas core is not 
in the north region and the massive subhalo is located in 
the south region. Some discrepancy between the results of 
simulations and observations would be caused by a choice 
of their initial conditions. For instance, they assumed 
that the gas temperature is almost isothermal, while re¬ 
cent studies (Okabe et al. 2014b) pointed out that the 
gas temperature universally drops off beyond ^ 0 . 5 r 2 oo- 
The intrinsic scatter of the gas density deviated from the 
universal profile (Okabe et al. 2014b) might be also im¬ 
portant. The observed configuration of the radio relics 
and X-ray surface brightness distributions are realized by 
other numerical simulations (Token et al. 1995; Roettiger 
et al. 1999). A head-on merger with a mass ratio of 1:4 
in less than IGyr after core passage (Token et al. 1995) 
have shown that the largest gas core and the prominent 
shock are found in the more and less massive halos, re¬ 
spectively. Roettiger et al. (1999) have conducted the 
magneto-hydrodynamics (MHD) and dark matter simu¬ 
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lation to study A3667 and shown that a cluster merger of 
a mass ratio of 1:5 reproduce a largest radio relic associ¬ 
ated with a subhalo and X-ray emission associated with a 
main halo. The mass ratios in two simulations are some¬ 
what larger than that for the cluster. Those simulatious 
assumed the first order Fermi acceleratiou to compute the 
radio relics. In order to precisely realize the observed fea¬ 
tures by numerical simulations, further studies using ob¬ 
servational results would be important to understand the 
formation of radio relics. Cosmological MHD simulations 
(Skillman et al. 2011, 2013) would be also important to 
understand the statistical properties of the radio relics. 

In contrast to giant radio relics, no prominent radio 
halo is found in CIZAJ2242.8-1-5301 (van Weeren et al. 
2010; Stroe et al. 2013) in the central region even at lower 
frequency 152MHz. That means that relativistic electrons 
are neither fully accelerated nor re-accelerated. For exam¬ 
ple, dying relativistic electrons escape from the backward 
side of the moving direction of radio relics, but are not 
re-accelerated by downstream turbulence. The fact leads 
to the condition 


Gcc ^ ^merger}; 


( 1 ) 


where face ^cooi and ^merger are the acceleration timescale, 
the radiative cooling timescale and the merger time 
scale related to the particle acceleration, respectively. 
The Torentz factor of relativistic electrons is 7 ~ 6.6 x 

10^ ( i 53 MHza+L 92 ) ) ' (ifj) ^ ■ Given the Torentz 
factor, the cooling time of relativistic electrons is com¬ 
puted by the synchrotron and the inverse Gompton losses, 
as follows. 


^c°o'-0-2(6xV) 

{ ^ V ( -SqMB \ ^ 
^ U-25(i + -2)2AiGy 


Gyr, ( 2 ) 


where is the energy density of the cosmic 

microwave background (GMB). We define the merger 
timescale, tmerger, as the time that the northern subhalo 
travels from the center-of-mass. 


t 


merger 


- 0.2 


( ^NM-CM ^ 

/ U N 

\440 /i“^kpc / 

V 2000 kms-iy 


Gyr, (3) 


where we assume that the merging velocity v = 
2000 km s“^ is equivalent to the shock velocity estimate 
by X-ray observation. We assume stochastic acceleration 
due to gyroresonant interaction with Alfven waves (e.g. 
Ohno et al. 2002; Fujita et al. 2003; O’Sullivan et al. 2009), 
in the framework of the second-order Fermi acceleration. 
The acceleration timescale is computed by 


tare — 


K, 


(4) 


where va is the Alfven velocity and Kg — ccL/^j-es is a 
spatial diffusion coefficient parallel to the mean magnetic 
field by a pitch angle scattering, respectively. Here, cl is 
the Tarmor radius and ^res = (^J?res/J?)^ is a fractional 










Radio Relics Tracing the Projected Mass Distribution in CIZA J2242.8+5301 


7 


No. ] 


energy of MHD turbulence at the wavelength satisfying a 
resonance scattering condition. For simplicity, we ignore 
the wave damping effect. The balance between the cooling 
time/merger time and the acceleration time constrains the 
fractional energy of the turbulence of 


< 8 X10 


-7 


B 

TJIg) 


-4 


B y 

TJIg) 


Bcmb 


3.25(1+ z)2/iG 


and 


e™-- < 6 X 10- 


Wg 


-7/2 


-1 


t^merger 


0.2Gyry 


We here assume the electron density Ue = 10“^ [cm“^] 
and the mean molecular weight /r = 0.59. The two con¬ 
straints of Cres°^ Comparable to each other. 

We note that the two equations above have different de¬ 
pendencies on the magnetic field strength, B. These con¬ 
straints change slightly depending on the choice of the B- 
field strength. Although the exact limits would depend on 
the spatial position, our observations provide a crude up¬ 
per limit on the fractional energy of the MHD turbulence 
as ^res < 0(10“®). It indicates that the MHD turbulence 
is significantly suppressed in the central region, at this 
merging stage of the cluster. In the above discussion, we 
ignore the mechanism to generate the Alfven waves. If 
the Alfven waves is connected by the cascade of the fluid 
turbulence, the balance between the energy loss and gain 
(Gcc = tcooi) does not change but Gcc = ^merger is modi¬ 
fied to face + tcas = ^merger, where teas is the time scale of 
the cascade. As a result, the upper limit on be¬ 
comes larger than that on We also note that the 

estimation of the MHD turbulence depends on the pre¬ 
cise mechanisms of particle acceleration on the basis of 
microplasma physics, and thus further studies regarding 
other various mechanisms are essential for a deeper under¬ 
standing of radio halos. If future deep radio observations 
discover very faint diffuse emission, the constrains on the 
MHD turbulence will be much improved. 

We compare our weak-lensing masses with previous 
study (Jee et al. 2015). Jee et al. (2015) performed a 
weak-lensing mass estimate for the cluster using a color- 
magnitude selected galaxy sample that contains cluster 
red-sequence galaxies fainter than 22 ABmag. However, it 
has been well established by earlier work (e.g. Broadhurst 
et al. 2005; Okabe et al. 2013; Okabe & Smith 2015) that 
such a background selection can lead to systematic dilu¬ 
tion of the weak-lensing signal, resulting in a substantial 
underestimate of the cluster mass especially at inner radii. 
By fixing the concentration parameter according to the 
mean c-M relation of Duffy et al. (2008) that is based on 
the WMAP five-year cosmology, they estimated M^qq = 
X 1014 H-^Mq and = 7.7^11 x lO^^ H-^Mq 

for the southern and northern halos, respectively. Their 
mass estimates for the northern and southern compo¬ 
nents, as well as their total mass estimate, are compat¬ 
ible with our results within large errors inherent in noisy 


weak-lensing measurements with a small number of back¬ 
ground galaxies. On the other hand, Jee et al. (2015) 
concluded that the cluster is the result of an equal-mass 
merger, whereas our best-fit model supports a ^ 2 : 1 ra¬ 
tio merger. In order to understand the origin of this dis¬ 
crepancy, we re-performed two-dimensional shear fitting 
assuming the single scaling concentration of Duffy et al. 
(2008) as done by Jee et al. (2015). We obtained MIqq = 
8.i;^;® X 1014 H-^Mq and x 10i4 H-^Mq 

for the southern and northern halos, respectively. The 
resulting mass ratio is thus reduced to ~ 1.2, in better 
agreement with their results. Therefore, their conclu¬ 
sions may be sensitive to the particular choice of the mean 
concentration-mass relation. 

7. Summary 

We have conducted weak-lensing analysis using three- 
band imaging of Subaru/Suprime-Cam. We have carefully 
selected background galaxy populations free of contami¬ 
nation by unlensed cluster galaxies, following the color- 
color selection method developed by (Medezinski et al. 

2010) . Two dimensional shear fitting enabled us to mea¬ 
sure masses for two halos and their positions (Okabe et al. 

2011) . We have found that the southern halo is the main 

cluster and the northern halo is the colliding subcluster 
of which mass is about half of the main cluster mass. 
Spatial offsets between luminous galaxies and the weak- 
lensing-determined centers are at maximum 2cr level. The 
mass-to-light ratios for two halos within from 

weak-lensing-determined centers are comparable because 
the concentration parameter for the subcluster is higher 
than that for the main cluster. 

The distance ratio between outer edges of the two relics 
and the center of mass are comparable to those between 
two halo positions and the center of mass, as well as be¬ 
tween the Mach numbers determined by X-ray observa¬ 
tions. The gas halo is associated with the southern main 
halo, while no well-defined gas core is found in the north¬ 
ern subhalo. Our mass measurement could consistently 
explain the merger configuration observed by X-ray, radio 
and weak-lensing techniques. The result implies a possi¬ 
bility that the radio relic is formed after the core passage. 

We have found that the curved shapes of radio relic 
coincides well with the projected mass distribution, in¬ 
dicating that radio relics on cluster-scale well follow the 
global structure of the mass distribution. 

No prominent radio halo is found, which enabled us to 
put an upper limit of the fractional energy of MHD turbu¬ 
lence at resonant wavenumber of {SB/B)"^ < 0(10“®), by 
comparing the acceleration timescale of the resonant ac¬ 
celeration with the cooling time or the merger time after 
the core passage. 
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